Probing the Dynamical Behaviour of Surface Dipoles Through Energy Absorption 

Interferometry 
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Spatial interferometry, based on the measurement of total absorbed power, can be used to de- 
termine the state of coherence of the electromagnetic field to which any energy-absorbing structure 
is sensitive. The measured coherence tensor can be diagonalized to give the amplitude, phase, po- 
larization patterns, and responsivities of the individual electromagnetic modes through which the 
structure can absorb energy. Because the electromagnetic modes are intimately related to dynamical 
modes of the system, information about collective excitations can be found. We present simulations, 
based on the Discrete Dipole Approximation (DDA), showing how the dynamical modes of systems 
of surface dipoles can be recovered. Interactions are taken into consideration, leading to long-range 
coherent phenomena, which are revealed by the method. The use of DDA enables the interferomet- 
ric response of a wide variety of objects to be modeled, from patterned photonic films to biological 
macromolecules. 

PACS numbers: 42. 25. Kb, 42. 25. Hz, 42.50.Ar, 44.40.-|-a 



I. INTRODUCTION 

Understanding the processes by which molecules and 
particles on surfaces absorb and scatter electromagnetic 
radiation is an area of considerable scientific importance. 
Information may be needed because of an intrinsic de- 
sire to study surface physics, because an active surface 
forms part of an electronic device, or because of the need 
to locate defects and identify contaminants. Examples 
include the dynamics of adsorbed biomolecules the 
adhesion and movement of droplets [2| , the plasmonic be- 
haviour of nanoparticles [1, |3| , and the photonic response 
of patterned films 0, @. Surface processes have tradi- 
tionally been studied through infrared and optical mi- 
croscopy, spectroscopy, and scattering, with AFM play- 
ing a key role in recent years. 

Often, the information needed has to be acquired in- 
directly through measurements made at inappropriate 
wavelengths, or through models that are hard to extrap- 
olate. Crucially, it is necessary to measure the state of 
coherence of the field to which a system is sensitive in or- 
der to fully characterize its behaviour: for example when 
calculating radiative heat transfer between nanoparticles 
0, II] . Measurements made with a single source only gain 
access to a subset of the information available. Of par- 
ticular importance are those cases where the ability of a 
surface to absorb energy is central to the operation of a 
device. For example, complex schemes are being explored 
to increase the efficiency of photovoltaics d, [13] , to op- 
timize the sensitivity of infrared and optical sensors , 
and to develop frequency selective surfaces, high-power 
mirrors, and hyperspectral black absorbers. In all cases, 
the physics of the absorption process is intimately related 
to the existence of coherent phenomena such as resonant 
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dipole interactions, surface plasmon coupling [l2!|, and 
surface phonon exchange. Even in the case of apparently 
simple structures, such as planar dielectric and metallic 
interfaces, the near-field spatial and polarimetric corre- 
lations contain a wealth of information p^l - fTsj . 

A typical set of questions is as follows: What are the 
scale sizes and forms of the dynamical modes responsi- 
ble for absorbing energy? How many modes are involved 
in the absorption process? What are their responsivi- 
ties? What physical mechanisms mediate the collective 
behaviour? In this paper we argue that it is possible 
to measure the state of coherence of the electromag- 
netic field to which any system is sensitive, and thereby 
determine the spatial forms, polarizations, and relative 
strengths of the individual modes responsible for absorb- 
ing energy. Because the electromagnetic modes are inti- 
mately related to dynamical modes of the system, it is 
possible to accrue direct information about the dynami- 
cal modes themselves, and thereby reveal the existence of 
collective phenomena. The basic technique can be used 
at any wavelength, and can be implemented in a number 
of different ways. 

The basic idea is to illuminate the Structure Under 
Test (SUT) with a pair of phase-locked, near-field or far- 
field point sources. As the differential phase between the 
sources is rotated, the total absorbed power displays a 
fringe. If the complex visibility of this fringe is mea- 
sured for different pairs of source locations, the dynami- 
cal forms of the modes responsible for absorbing the en- 
ergy can be found. In the case of detectors and pho- 
tovoltaics, the scheme is particularly easy to implement 
because the degrees of freedom responsible for generating 
the measured output are those of practical significance. 

Although the technique is applicable to a wide range 
of problems, we illustrate the method by showing how 
it provides information about the collective behaviour of 
systems of surface dipoles. The SUT is modeled by using 
the Discrete Dipole Approximation (DDA) [16l - [l8| . to- 
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gether with a polarizability that aUows for bound states 
and damping to the substrate. DDA was originahy de- 
vised [l^ to model hght extinction by astrophysical par- 
ticles, but has since been developed in areas as diverse as 
mode ling the aggregation of irregular particles from col- 
loid s [2011 , the scattering of light by particles on surfaces 
[21I |22|. the scattering of light by periodic targets [l^ , 
and the optical properties of viruses, proteins, DNA, and 
blood cells |24l426l |. Formulating our measurement pro- 
cedure in terms of DDA opens the way to predicting the 
interferometric response of a wide variety of complex ob- 
jects such as biomolecules, droplets, and nanoparticles, 
and thereby creating representative dynamical models of 
measured interferometric data. Other simulation meth- 
ods, such as MoM, can be used. 

As an aside, a variant of EAI can be used with DDA 
purely as a numerical tool for determining the state of 
coherence, and therefore the incident electromagnetic 
modes, to which any irregular dielectric body is sensitive. 
The scheme differs from the usual method for determin- 
ing the state of coherence of thermally radiated fields, 
through 



W{r,r') = /iofe^^coth 



ImG (r,r'), (1) 



where G (r, r') is the electric Green's dyadic of the whole 
system including scattering, because ((!]) requires the ob- 
ject to be in thermal equilibrium with its surroundings, 
and therefore includes the blackbody environment and 
scattered fields. As a consequence, even lossless dielectric 
objects display near-field coherence efi'ects, as discussed 
by Yannopapas in the context of chains of spheres [27| . 
In contrast, EAI only includes those degrees of freedom 
responsible for absorbing energy. 

In section |ll] we present the theory of Energy Absorp- 
tion Interferometry (EAI), and then in section |lll] illus- 
trate how it can be used to recover the dynamical modes 
of ring-like and linear chains of coupled dipoles, with and 
without defects. Linear chains of dipoles (28j|, spheres 
[27j . and cylinders [2^ are of practical importance, and 
exhibit complex forms of spatial and spectral coherence. 
We show how EAI can reveal large-scale collective be- 
haviour, and how the response of a discrete system tran- 
sitions to that of a continuous system when the SUT 
contains many dipoles per wavelength. In section HVl we 
discuss practical implementation. 



II. BASIC SCHEME 

A. Dynamical Behaviour 

Consider a collection of particles whose electromag- 
netic behaviour can be modeled as a system of dipoles. 
Although EAI can be applied to any three-dimensional 
distribution, assume that the particles are bound to a 
surface. Fig. [TJ The dynamical response of dipole i can 
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FIG. 1. A collection of particles on a surface. The individual 
particles have dipole moments pi. The system is illuminated 
by a pair of phase-locked sources, which can have any polar- 
isation and position. The experiment measures the complex 
visibility of the absorbed power as the differential phase be- 
tween the sources is varied. 



be described by the Lorentz model 



at 



-UJlMt)^—'Eir^,t)■i^i^, (2) 

TO,: 



where is the effective separation vector of charge qi 
from its partner pinned at r.^. The dyadic product x^Xi 
constrains qi to move in the Xi direction, but other dyads 
could be used to place more general constraints on the de- 
grees of freedom of motion, ujoi is the resonant frequency 
due to binding, rrii the mass, and F.^ the momentum re- 
laxation rate of dipole i. The damping arises from a 
transfer of energy to the phonon or electron system of the 
substrate. ^ does not include a radiation-reaction term, 
because we are only interested in the absorbed power, 
and the radiation-reaction term has little effect on be- 
haviour. E(ri,t) is the applied field, which includes the 
incident field, E^(ri,i), and the scattered field from all 
other dipoles, E'^(ri,i). In the parlance of DDA, E(rj,t) 
is the 'local field' because it does not include the self field 
of the dipole at r^. Invoking time-harmonic behaviour. 



, leads to a polarization of 



P(r) 



V "^^ ('^o, 



-E(r,) • XiXj 



(3) 



The anisotropic susceptibility tensor, Xj(r), defined 
through P(r) = eoXi(r) • E(r), becomes 



(4) 



It will be convenient to describe behaviour in terms of 
polarisation current density J(r) = dP{t)/dt, giving 
J(r) = ^(r) ■ E(r), where 



^(r) = E 



qf -iuj5{r-ri) 
TO, (w^j - a;2 - iuiTi)' 



(5) 



Large currents can flow at the resonance frequencies, ujoi, 
and an increase in either binding, tOoi, or damping, F^, 



3 



decreases the low-frequency conductivity. The permittiv- 
ity and conductivity tensors are vaUd in both the high- 
frequency and quasi-static regimes. 

To simulate EAI, it is necessary to calculate the total 
time-averaged power, W, transferred to the substrate. 
The power can be calculated by treating the system as a 
dielectric or conductor Adopting the latter, 



W = ^Re J E*(r) ■3{r)d^r 



(6) 



Re / E*(r) • ^(r) • E(r) d^r, 



giving 



2m, 



1 



|E(r.)-x,|^ (7) 



in the symmetric form 
Re I / E*(ri) •^(ri,r2) •E(r2)d'rid3r2, (8) 



It is convenient to present 
where 



o-(ri,r2) = cr(ri)^(ri - rz) 



(9) 



is diagonal in position, which emphasizes its local form. 
When applied to the system of dipoles, 

W^lRcJ2J2^*ir^)■^^^'^i^J)^ (10) 



where 



if 



—iuj 5j 



rrii {loI 



iujTiY 



(11) 



is a discretized conductivity tensor. 

([6]) and ([8]) cannot be used unless the total electric 
field, E(r) is known, but the total electric field is the 
sum of the illuminating field, E^(r), and the scattered 
field, E'^(r), giving 



J(r)=^(r). [E^(r)-HE^(r)] 



(12) 



The scattered field can be calculated through the Green's 
dyadic of the environment. 



E^(r) = / G(r;r') • J(r')d^r' 



(13) 



where all multiplicative factors have been included in 
G(r;r'). Green's dyadics of dipoles on, and under, di- 
electric surfaces have been studied extensively, including 
multilayers [sO] and surface roughness [Slj. Ordinarily, 
we might use the Green's dyadic corresponding to the 
dielectric half space shown in Fig. [1] but for illustrative 
purposes, and consistency with DDA, we shall use the 
Green's dyadic of free space, ignoring the complex re- 
fractive index of the substrate. Of course, the surface 



could also be modelled through DDA. In the case of par- 
ticles on very thin, say SiN, membranes, this assumption 
is valid without further work, but in any case, it does 
not alter the substance of what follows. The same func- 
tional forms shall be used to describe scattering by the 
dipoles and the excitation of the system by external point 
sources. The free-space Green's dyadic is 



G(r;r') = 



.9(r|r'), 



where g{Y\ r') is the scalar Green's function 
1 



g{r;r') = 



47r r — r' 



-exp [ifco|r-r'| 



(14) 



(15) 



V denotes the gradient with respect to the observation 
point r. The first term in (jl4p emerges from the magnetic 
vector potential, and the second term from the electric 

scalar potential, and therefore G(r; r') includes coupling 
due to current flow and the exposure of charge. 

In order to perform numerical simulations, it is conve- 
nient to express G(r;r') in Cartesian coordinates. The 
Green's dyadic has near-field (NF) i? < A, intermediate- 
field (IF) i? ~ A, and far-field (FF) i? > A, contributions: 



G 



'NF 



G 



IF 



'FF- 



Setting R — r j — , and evaluating (fT4| gives 
exp [ikR] 1 



Gnf 

GiF 



exp [ikR] i 



= 3RR 



i?2 



(16) 



(17) 



GpF = iijJiiQ 



inR kR 
exp [ikR] 
AttR 



3RR 



i?2 



RR 



where R = |R|. By expressing the Green's dyadic in this 
way, individual terms can be turned on and off during 
simulations to gain a better understanding of behaviour. 

It can be shown that as w ^ 0, Gatf leads to the quasi- 
static polarization of the dipoles, and therefore electro- 
static coupling is included. We do not include the depo- 
larization term —i{uifj.o/3k^)lS(R), consistent with the 
use of the localized rather than total field. In DDA, 
where each dipole represents a small cell of a large di- 
electric body, care in needed to distinguish between the 
total and localized fields [l3| , particularly if power is dis- 
sipated in the dielectric itself. EAI is equally applicable 
to this case, but for presentational purposes, we have as- 
sumed that each particle is small enough that it can be 
represented by a discrete dipole. 

There are several possible ways forward: substitute 
(IT51) in (|12l) and attempt to solve for J(r), or substitute 
(IT2t in (IT3l) and attempt to solve for E(r). The second 
approach is well conditioned, and so 



E(r) = E^(r) +E^(r) 



(18) 



E^(r) + / G(r; r') • ^(r') • E(r') rf^r'. 
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Then 



where, symbohcally, 
corresponding to 



E(r') = / M(r;r') •E^(r)d='r 



(19) 



M(r; r') 



is the inverse kernel 



M(r;r')= K(r - r') - G(r; r') • ^(r') 



(20) 



Because only a finite number of dipoles is involved, the 
inverse operator ([T9|) exists. Substituting CT(r') in ([T8|. 
and noting that it is only necessary to calculate the total 
field at discrete points in order to calculate the dissipated 
power, 



^M(r,;r,)-E(r,) = E^(r,), 



(21) 



where 

r ? 

M(r,;r.) = ^K, - _ _ .^r.) ^^^^ ' '''^ ' 

(22)' 

In DDA, parameters in the Green's and conductivity ten- 
sors are combined to express the scattered field in term 
of a modified Green's tensor and the susceptibility ten- 
sor. The physics is, however, the same. It is convenient 
to write (l22]l in the form 



1 - UJi 



-G (rf,ri) ■ XiXi 



(23) 

where the tilde denotes normalization with respect to the 
resonant frequency of the bound dipole: uji = uj/uJoi and 

Ti = Ti/ujoi- G {rj^Ti) is the Green's dyadic introduced 
previously, but with the prefactor uj^q changed to k^. 
According to (jlj, = qf/{miujl^) is the electrostatic, 
w — 0, polarizability of dipole i. The Clausius-Mossotti 
relationship can be used for simple shapes: 



i^r - 1) 

l + n{er - 1) 



V, 



(24) 



where V is the volume of the particle, having shape factor 
< n < 1. 

Because the scheme is fully discretized, it is convenient 
to introduce matrix notation. If and e are column 
vectors containing the complex amplitudes of the incident 
and total Cartesian field components at the positions of 

A-^e'. In 



(25) 



the dipoles, then (|21l) can be solved: e 
matrix notation, (jlOp becomes 



W : 



iRe[etZe] 
2 



= - e^t U 
2 L 

= e^tLe^ 



-11 1 



where Y. is the matrix equivalent of (|TT|) . and 51'' is the 
real part. Taking the trace of each side of the last line of 
(|25p . noting that the trace is invariant to cyclic permu- 
tation, and calculating the expectation value gives 



(W^) =Tr[(e^e^t)L] = Tr [E^L] . 



(26) 



is Hermitian, and contains the correlations between 
the Cartesian components of the incident field at the lo- 
cations of the dipoles. gives the expectation value 
of the absorbed power in terms of the state of coherence 
of the incident field and a response matrix L. Not only 
is ()26p computationally convenient, it is conceptually im- 
portant. The trace of the product of a matrix with the 
adjoint of another is an inner product in the abstract vec- 
tor space of complex matrices. (pS)) therefore describes 
the projection of the state of coherence of the incident 
field onto the state of coherence of the field to which the 
system is responsive. In the case of continuous absorbers, 
an integral formulation of the same concept is available 
[sll. E^ and L are Hermitian, and can be diagonalized: 
= Y^m^mam^ln, ^ ^ J2n ^nPn^l, which Can bc Sub- 
stituted in (l26l) to reveal 



(27) 



(l?fl) describes the way in which the modes of the incident 
field v,„ project onto the modes of the system, u„. The 
am are proportional to the levels of excitation of each 
of the incident modes, whereas the /3„ correspond to the 
relative responsivities with which the system absorbs en- 
ergy in each of its modes. According to (l^5|) . the response 
matrix L contains two effects: the intrinsic ability of the 
dipoles to absorb energy, Z"", and the appearance of col- 
lective behaviour due to scattering, [M]~^. 

The electromagnetic Green's dyadic induces a wide va- 
riety of coherent behaviour, but additionally an elastic 
Green's dyadic could have been included in ^ to de- 
scribe acoustic coupling through the excitation of sur- 
face waves. The dynamical modes of interest might then 
be mediated by long-range phonon exchange. The over- 
all principle, however, remains the same. Also, because 
DDA has proven to be an effective tool for modelling 
electromagnetic scattering from complex objects, includ- 
ing internal retardation, it use here for describing EAI 
means that the interferometric response of macroscopic 
inhomogenious objects, such as particulate material or 
droplets of a liquid, can be calculated. In this case, each 
macroscopic object is represented by a large collection of 
dipoles, rather than only having one dipole per object. 



B. Interferometry 

In an experiment, the SUT is illuminated by a pair of 
phase-locked sources, and the absorbed power recorded 
as the differential phase between the sources is rotated. 
The absorbed power displays a fringe, and the complex 
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visibility of this fringe is measured. In reality, the fre- 
quencies of the sources can be offset slightly, and the am- 
plitude and phase of the fringe measured directly by using 
lock-in techniques. If the complex visibility is recorded 
for different pairs of source locations, either in the near 
field or far field, or a mixture of both, the response matrix 
L can be found. 

Suppose that two sources are used, one producing the 
incident field e{e-''^^ , and the other producing ejC-''^^ . Ac- 
cording to (|25p , when both sources illuminate the sample 
simultaneously. 



W ■ 



■ft - 



L[e(< 



(28) 



Ensemble averaging is not required because the sources 
are fully coherent. ([25)1 has four terms. 



= 1^11 1^22 + 1^126' 



(29) 



W2ie 



where A0 = 02 ^ 0i is the phase difference between the 
sources. The first two terms, Wu and W22, are real, and 
give the powers dissipated by the sources individually. 
The second two terms are complex conjugates of each 
other 14^12 — W21 , and produce a fringe when the differen- 
tial phase is rotated. It follows that W^n, W22, W^2i, W^i2 
can be found experimentally. If the field vectors efe^'^^ 
and 626-''^^ corresponded to point sources at the positions 
of the dipoles, the response matrix L could be found di- 
rectly, but in reality only distributed fields are possible. 

Although the response matrix can be computed for 
different forms of illumination, plane waves, spherical 
waves, Gaussian beams, etc., here we assume that point 
sources, in the form of miniature current probes, are used. 
The current density of a point source at is 



J^(r') =(5(r'-r^)I^, 

where I'^ is the vector- valued current-moment, 
source produces an incident field of the form 



E^(r) = 



G''(r;r') 



J^(r')d'r', 



(30) 
This 

(31) 



The dissipated power can be written 



W= iRe[eti^] 



(34) 



^tG^t [M-i]^Z'-[M-i] G^i 



P-.P 



L is the intrinsic response matrix with respect to the in- 
cident field at the positions of the dipoles, whereas H is 
the response matrix with respect to current moments at 
the positions of the external sources. 

To carry out interferometry, two sources are needed. 
Imagine that the current moment corresponding to 
source 1 is \ fe^'^^ , and the current moment corresponding 
to source 2 is if e^''^^ . if and if are unit vectors having a 
single real entry corresponding to a single current source 
with a single Cartesian component. When both sources 
illuminate the sample simultaneously. 



W : 



•-Pf - 



I2 'e 



H 



fife^ 



+ h e 



Pf,j4>2 



(35) 



Ensemble averaging is not required because the sources 
are coherent. (j35p has four terms 



W = if^Hi( -r ,2 

= Hii + H22 



^tHif + iftHif + ifTHi( 



Hi2e- 



jA0 



H2ie-^'^'^. 



(36) 



The first two terms, Hn and i?22, are real, and give 
the powers dissipated by the sources individually. The 
second two terms satisfy H12 = -ffliJ and produce a 
fringe when the differential phase is rotated. Thus by 
measuring the power absorbed when each source is ap- 
plied separately, and the amplitude and phase of the 
fringe when the two sources are applied simultaneously, 
Hu, H22, H21, and H12 can be found. 

Now move the two sources around within some sam- 
pling volume, stepping from one sample position to an- 
other. Usually the sample points will lie on a line or plane 
in the near field or far field, but more complicated paths 
can be used depending on the sampling needed. Con- 
struct the matrix 1^, where column n, if, corresponds to 
the source being at position n. Guided by form the 
new matrix 



where G (r; r') is the Green's dyadic that couples the 
source current to the incident field at the dipoles. Sub- 
stituting (|30l) in (|3T1) leads to the discretized form 



(32) 



For a single, simply-polarized source, i^ has one non- 
zero entry corresponding to the position and polarization 
of the probe. Using ([32]) and calculating the induced 
current moment, 
gives 



at each dipole by integrating 



-iG^i^. 



(33) 



Q-l^tG^^t [M-i]Tz[M-i] G^l 



PiP 



(37) 



Q is not Hermitian, H ^ Y.^ . Consider the Hermitian 
part of Q, 

i [Q + Qt] = l^tG^'t [M-i] ^ r [M-i] G^l^ (38) 
= l^tG^'tLG^iP 

According to ([38]), the n'th element of the leading diag- 
onal of H corresponds to the total power absorbed when 
single source n illuminates the sample. The off-diagonal 
terms describe the fringe when a pair of phase coherent 
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sources n and n' are used. It is therefore possible to pop- 
ulate H experimentally. If one source is held fixed at n, 
and the other moved through all other positions n', the 
complex visibility of the spatial fringe is 



2H„ 



In 



(39) 



It is possible to plot out |7„„'|, with respect to reference 
point n, by moving along the n'th row of H. These plots 
reveal the transverse and longitudinal coherence lengths, 
areas, and volumes of the polarized point-source excita- 
tion to which the structure is sensitive. H is the response 
matrix of the system with respect to point source excita- 
tion in the scanning volume. Once it has been populated 
experimentally, it can be diagonalized: 



H = G^tLG^ 



iln 



(40) 



1.5 



o ' 

Q. 

■o 

O „ ^ 

« 0.5 





1 1 

ll 


1 1 


- 


'1 ; 
l| - 
ll :: 






'1 
1 1 






1 1 

' \ 

- — 1^ • - 


A' / " 



100 



200 300 
Frequency (GHz) 



400 



500 



FIG. 2. Total power absorbed (arbitrary units) as a function 
of frequency for two parallel dipoles, 0.1 mm apart, each hav- 
ing ujo = 300 GHz, and F = 20 GHz. The lines correspond to 
polarizabilities of (solid) a/eo = 0.0005 mm^, (dotted) a/eo 
= 0.005 mm^ (dashed) a/to = 0.01 mm^ 



Substituting PO)) in the last line of ([M)) . shows that if 
the source vector corresponds to each of the normalized 
eigenvectors in turn, = w^, then the eigenvalues are 
the relative responsivities of the system. The corre- 
spond to source vectors that excite modes of the system. 
Following ((33)) . the dipole moment vector of mode m 
is 



1 



(41) 



No attempt is made to recover the displacements, as gen- 
erally the charge at each site is not known. 

Response matrix H combines three effects: (i) The gen- 
eration of distributed fields by point current sources, de- 
scribed by G'^; (ii) scattering interactions between the 
dipoles, described by [M~^]; and (iii) the intrinsic abil- 
ity of each dipole to absorb energy, described by T.^. The 
second and third contributions, which generally cannot 
be separated, are intrinsic to the system, whereas the 
first term is a feature of the measurement. One could 
calculate 



(42) 



and although these are field modes at the positions of the 
dipoles, they still relate to the original source volume. 
To calculate the modes at the positions of the dipoles, it 
is necessary to 'deconvolve' the illumination patterns of 
the sources. Using (PH)) . the deconvolved field modes are 
found by diagonalising 



(43) 



Generally, the Green's matrix will be singular, or poorly 
conditioned, depending on how well the sample points 
and polarizations span the vector space of fields at the 
dipoles. For poorly conditioned mappings, the general- 
ized inverse should be used to recover those degrees of 



freedom that are available to the experiment. In reality, 
it is only necessary to include enough sample points to 
determine the significant degrees of freedom contained in 
L, which may be fewer than the number of dipoles. One 
approach is to keep on adding, during an experiment, 
more and more points until the measurement has con- 
verged. Incremental Singular Value Decomposition can 
be used to diagonalise L in real time until all of the avail- 
able degrees of freedom have been found. Finally, the 
incident fields can be mapped onto the dynamical modes 
of the dipoles through 



1 



(44) 



III. SIMULATIONS 

A typical simulation of an experiment proceeds as fol- 
lows. Establish a system of dipoles, and calculate the 
conductivity matrix Y. using (fTTj). Calculate the scatter- 
ing operator M using ([23l) . and invert it. Calculate L as 
defined by and diagonalize it to find the eigenval- 

ues and eigenvectors u„ and /?„. These are the collective 
dynamical modes of the system responsible for absorbing 
energy, which we would like to find. Choose a selection 
of source positions and polarizations, calculate the source 
Green's matrix G^, and the response matrix H through 
(HUl) . Determine the amplitudes and phases that would 
be measured experimentally using p4l) . reverse the pro- 
cess using (j43|) . and calculate the apparent dynamical 
modes using (l44t . Compare the actual modes with the 
recovered modes. Information may be lost due to the 
sampling used, the introduction of experimental errors, 
and the complete loss of certain fringe visibilities due to 
noise. During the work we found it instructive to produce 
'movies' showing the time evolution of the actual and re- 
covered modes. It should be stressed that in all of the 
simulations that follow, only the total power absorbed by 
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FIG. 3. Total power absorbed (arbitrary units) as a single 
source is swept in the a;-direction, at a distance y — 5 mm from 
two parallel dipoles, 0.1 mm apart, each having cjo = 300 GHz, 
r = 20 GHz, and a/eo ~ 0.005 mm^. The source frequen- 
cies were (dashed) 240 GHz, (solid) 300 GHz, and (dotted) 
340 GHz. 



the complete system was regarded as measurable. 

First, consider the resonant coupling of two parallel 
dipoles. Figured] shows the total power absorbed, in arbi- 
trary units, as a function of frequency for two z-polarized 
dipoles, placed 0.1 mm apart on the x-axis, each having 
ujo — 300 GHz, and F = 20 GHz. The z-polarized source 
was located at x—y—10 mm to ensure that symmetric 
and antisymmetric modes could be excited. The lines 
correspond to different polarizabilities, thereby changing 
the degree of coupling, and the individual conductivities. 
The solid line shows the case when the polarizability is 
low, a/eo = 0.0005 mm"^. The dipoles respond indepen- 
dently, giving rise to a single spectral feature at 300 GHz, 
having a width of F = 20 GHz. Because the conductiv- 
ity is low, the absorbed power is small. The dotted line 
shows the case when the polarizability is increased to 
a/eo = 0.005 mm'^. The dipoles interact, and split into 
an antisymmetric mode at 240 GHz, having a width of 
r/2, and a symmetric mode at 340 GHz, having a width 
of 2F. The delocalized excitation is mediated mostly by 
near-field coupling, and the splitting is caused by a de- 
crease (increase) in the effective spring constant of the an- 
tisymmetric (symmetric) mode as a consequence of the 
interaction. The absorbed power is larger than before 
because of the increased conductivity. The dashed line 
shows the case when the polarizability is high, a/eo = 
0.01 mm'^, and the splitting large. The relative strength 
of the peaks depends on the position of the source, indi- 
cating that spectral information alone can be misleading. 
Although not shown, when the polarizability is increased 
to a/eo — 0.05 mm"^, the 'high-conductivity' limit is ap- 
proached, and the absorbed power falls to a very low 
level as the dipoles screen each other: the scattered field 
at each dipole is equal in magnitude and opposite in sign 
to the incident field. 

Figure |3] shows the total power absorbed as a single 
source (using the near-field Green's dyadic) is swept in 
the x-direction, at a distance y = 5 mm from the dipoles. 



with a/eo = 0.005 mm^. According to Fig. [2l the anti- 
symmetric mode peaks at 240 GHz, whereas the sym- 
metric mode peaks at 340GHz. The plots in Fig.[3]corre- 
spond to source frequencies of (dashed) 240 GHz, (solid) 
300 GHz, and (dotted) 340 GHz. Ahhough the dipoles 
are not resolved spatially at this scan distance, the anti- 
symmetric mode has a null on-axis, whereas the sym- 
metric mode has a peak on axis. From these absorption 
patterns one could infer the nature of the modes present 
at a single, intermediate frequency, say 300 GHz, which 
is the solid line in Fig. [3J but what happens if wish to 
determine the modal content directly using only a single- 
frequency source. To this end, two sources were held at 
fixed positions (source 1 at x = mm, y — 5 mm, and 
source 2 at a; = 5 mm, y = mm), and the differen- 
tial phase rotated. The fringe visibility was recovered 
from the simulation, and used to assemble the H matrix, 
which was then converted to the L matrix by 'deconvolv- 
ing' the beam patterns of the sources, and then diagonal- 
ized to find the natural modes and their responsivities. 
The simulated experiment was repeated for the polariz- 
abilities and frequencies used in Figs. [5] and |31 In all 
cases, the spatial dynamical forms and responsivities of 
the anti-symmetric and symmetric modes recovered were 
precisely the same as those calculated directly. The same 
results were achieved regardless of whether the source 
probes were described by the near-field, far-field, or full 
Green's dyadics. At 240 GHz and a/eo =0.005 mm^, 
the technique revealed a dominant anti-symmetric mode, 
corresponding to the peak in Fig. [51 At 300 GHz, the 
modes were present in roughly equal amounts, whereas 
at 340 GHz, the symmetric mode was dominant. Cru- 
cially, the forms and relative responsivities of the two 
modes present at 300 GHz could be determined using a 
simulated interferometric measurement, even though the 
local absorption spectrum does not indicate the presence 
of two modes. In the case where a/eo =0.00005 mm^, and 
over a wide range of frequencies, EAI correctly revealed 
the existence of two degenerate modes, corresponding to 
the oscillation of each dipole independently. Similar sim- 
ulations were carried out for different source types and 
locations, and the dynamical forms and responsivities of 
the two modes could always be recovered. This simple 
example illustrates the nature of EAI, but care is required 
because the method is based on the ability to 'deconvolve' 
the beam patterns of the sources, and so the accuracy of 
the final result will depend on the precision with which 
the experiment can be carried out. 

Next, a series of simulations were performed on ring- 
like configurations of parallel dipoles. By using rings, 
boundary conditions at the ends are eliminated, and dis- 
persion relationships can be found. Figure|3]shows the to- 
tal absorbed power as a function of frequency for the fol- 
lowing cases: (solid) 3 dipoles having a/eo = 0.005 mm"^ 
placed 0.1 mm apart at the corners of a triangle; (dotted) 
4 dipoles having a/eo = 0.003 mm'^ placed 0.1 mm apart 
at the corners of a square; (dashed) 8 dipoles having a/eo 
= 0.0004 mm^^ placed 0.0383 mm apart on the corners of 
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FIG. 4. Total power absorbed (arbitrary units) as a func- 
tion of frequency for a variety of ring-like configurations: 
(solid) 3 dipoles having a/eo = 0.005 mm'' placed 0.1 mm 
apart at the corners of a triangle; (dotted) 4 dipoles having 
a/eo = 0.003 mm^ placed 0.1 mm apart at the corners of a 
square; (dashed) 8 dipoles having a/eo = 0.0004 mm^ placed 
0.0383 mm apart at the corners of an octagon. The power for 
the 8-dipole case was scaled by a factor of 10 before plotting. 
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an octagon. The polarizabilities were chosen to give simi- 
lar spectral splittings, but the power for the 8-dipole case 
was scaled by a factor of 10 before plotting. The source 
was positioned at x =10 mm, y =0 mm, and operated in 
near-field mode. The absorbed power revealed a low-level 
modulation as the source was rotated around the center 
at fixed distance, but the polar patterns are not shown 
for brevity. In all cases, the spectra acquired 2 peaks as 
a consequence of resonant interaction. Interactions take 
place not only around the ring, but also across the ring. 

In the case of 3 dipoles, the peak at 378 GHz has a 
single dominant eigenvalue, corresponding to all dipoles 
oscillating in phase with equal amplitude. The peak at 
240 GHz, is two-fold degenerate and accounts for the 
symmetric and antisymmetric motion of the 2 remaining 
degrees of freedom. In the case of 4 dipoles, the largest 
eigenvalue of the peak at 357 GHz corresponds to the 
collective delocalized in-phase oscillation of all dipoles, 
whereas the peak at 289 GHz is 2-fold degenerate and ac- 
counts for opposite pairs being out of phase. The system 
of 8 dipoles produces 2 peaks, with the one at 454 GHz 
being the in-phase collective oscillation of the whole sys- 
tem. The one at 390 GHz is two-fold degenerate, and cor- 
responds to sinusoidal and cosinusoidal standing waves 
on the ring. 

Even by comparing the full spectra in Fig. 21 it is not 
possible to gain much information about the nature of the 
underlying modes. Moreover, the relative strengths of 
the peaks depends on the position of the source. If one is 
only interested in, and has access to, a limited frequency 
range, then there is little information in the polar power 
patterns that helps understand the number and forms of 
the modes. EAI can, however, reveal this information. 
To simulate EAI, two probes were rotated around the 
center of the system at a fixed distance of 5 mm, and 
the complex fringe visibility calculated for certain source 



FIG. 5. (a) Amplitude and (b) phase of the continuous fringe 
visibility of (solid) the triangular system at 240 GHz, (dotted) 
the square system at 289 GHz, and (dashed) the octagonal 
system at 390 GHz. The sources were swept around in a 
circle 5 mm from the centre. The reference angle was 90°. 



pairs. Fig. [5] shows the continuous forms of typical visi- 
bilities. Remember that these are not spatial fringe pat- 
terns, but the complex visibility as a function of angle 
for a particular reference angle. The various lines are 
for the (solid) triangular, (dotted) square, and (dashed) 
octagonal arrangements describe previously. The refer- 
ence angle was 90° in all cases, corresponding to a vertex 
point of the triangular and octagonal arrangements, and 
a face of the square arrangement. Figure [5] shows that the 
peak visibility is reduced when the sources are on oppo- 
site sides due to the 'leading' edge screening the 'trailing' 
edge. To recover the modes, a small number of sample 
points was used in each case; typically one more than the 
number of dipoles, avoiding the obvious degeneracies. It 
was found that the forms of the modes and their rela- 
tive responsivities could be recovered easily at randomly 
chosen frequencies, without much care in selecting the 
source positions. Also, the recovery process depended 
little on whether the sources were described by the near- 
field, far-field, or full Green's dyadics. Of course, only the 
full Green's dyadic describes the physical reality, but the 
tolerance to the precise form of the illumination pattern 
indicates that, in principle, it should be possible to carry 
out interferometry using a variety of different near-field 
and far-field probes. 

The octagonal system was considered again, but 
with the resonant frequency of one dipole increased to 
400 GHz, corresponding to a reduction in the mass of 
the oscillating charge. Although not shown here, the 
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FIG. 6. (Color online) Modal forms of the 413 GHz resonance 
of octagonal systems of dipoles with the resonant frequency 
of one dipole increased to 400 GHz. The reduced-mass dipole 
was moved around the chain (solid) none, (dashed) position 
1, (dotted) position 2, (dot-dashed) position 3. 



principle effect on the spectrum is to split the original 
low-frequency line at 390 GHz into two separate frequen- 
cies, 387 GHz and 413 GHz, reflecting the breaking of 
the degeneracy of the sinusoidal and cosinusoidal modes. 
There is also a quasi-static mode at 290 GHz. The domi- 
nant modal forms are shown in Fig. [51 at 413 GHz, as the 
reduced mass is moved around the octagon. In each case, 
the partner mode at 387 GHz has a standing wave that 
leaves the dipole with the reduced mass stationary. The 
corresponding visibility functions are shown in Fig. [7J 
where the reference angle was 90°. The visibility plots, 
particularly the phase, clearly show the movement of the 
defect dipole. 

A number of discrete source positions where then used 
in a attempt to detect the movement of the defect dipole 
around the octagonal structure. The sample points were 
placed on a circle, 5 mm from the center, 20° apart, with 
the first one opposite a vertex. The modes and their rel- 
ative responsivities could be recovered precisely regard- 
less of whether the near-field or far-field source Green's 
dyadic was used. The forms of most high-order modes 
could be recovered precisely with 8 or greater sample 
points, with some low-level of aliasing with 5 to 7 sample 
points, and not at all with 4 sample points. The determi- 
nation of the position of the defect dipole, depended little 
on the positions of the sample points. Overall, the modes 
could be recovered in detail despite the partial screening 
of the 'dark side' of the ring. 

Turn to linear arrays. Figure [5] shows the total power 
absorbed, in arbitrary units, as a function of frequency 
for 5, z-polarized dipoles, placed 0.1 mm apart on the x 
axis, each having lOq = 300 GHz, and F = 20 GHz. The z- 
polarized source was located at a:;=10 mm, y^lO mm to 
ensure that symmetric and antisymmetric modes could 
be excited. The solid line shows the case when the po- 
larizability is low, a/eg = 0.0005 mm^, and the dipoles 
are largely, but not completely, decoupled. The dotted 
line shows the case when the polarizability is increased 
to a/eo = 0.005 mm'^. The dipoles interact, and give 
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FIG. 7. (a) Amplitude and (b) phase of the recorded fringe 
visibility at 413 GHz of octagonal systems, with the resonant 
frequency of one dipole changed to 400 GHz. The lines cor- 
respond to (solid) none, (dashed) position 1 (corresponding 
to the 90° reference angle), (dotted) position 2, (dot-dashed) 
position 3. 
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FIG. 8. Total power absorbed (arbitrary units) as a function 
of frequency for a linear chain of 5 parallel dipoles, 0.1 mm 
apart, each having ujo ~ 300 GHz, and F = 20 GHz. The 
different lines correspond to polarizabilities, (solid) a/eo ~ 
0.0005 mm^, (dotted) a/eo = 0.005 mm , (dashed) a/eo = 
0.01 mm^. 



rise to 5 spectral features at 187, 239, 293, 332, and 
381 GHz; these alternate between symmetric and anti- 
symmetric forms. It is not necessarily the case that the 
number of features is equal to or smaller than the number 
of dipoles, because for large systems spatially harmonic 
response due to widely separated dipoles may exist. The 
dashed line shows the case when the polarizability is in- 
creased further, a/eo = 0.01 mm'^. The response becomes 
very similar to that of Fig. [21 but now both modes are an- 
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FIG. 9. (Color online) Modal forms of the resonances at 
(solid) 239 GHz, (dashed) 293 GHz, and (dotted) 332 GHz 
in Fig. [S] All dipoles are either perfectly in phase or out of 
phase. The lines are drawn to guide the eye only. 



tisymmetric whereas before one was symmetric and the 
other antisymmetric. Figure [S] shows the forms of the 
dominating modes, for a/eo = 0.005 mm'^, at the reso- 
nance frequencies of (solid) 239 GHz, (dashed) 293 GHz, 
and (dotted) 332 GHz. Each filled circle shows the posi- 
tion of the displaced charge at i = 0: Zi(0). Dipoles hav- 
ing the same sign in the figure remain perfectly in phase 
as time progresses, whereas those having opposite sign re- 
main perfectly out of phase. Those at zero remain at zero. 
The lines are drawn to guide the eye only. Table |T] gives 
the elements of the H matrix recovered from the fringes 
when the 5-dipole structure, centred on x =0.0 mm, was 
illuminated by 2 z-polarized sources at 280GHz, which 
according to the dotted line of Fig. [8] means that the 
sample was not resonant in any of its natural modes. 
A total of 5 source positions were used, x —\, 2, 3, 4, 
5 mm, with y = 5 mm in all cases. This set of data was 
then corrected for the beam patterns, and the dynami- 
cal modes calculated. The eigenvalues of the recovered 
modes are given in the next to last line of Table ID The 
forms of the recovered modes were then reconstructed 
and compared with the modes associated with the peaks 
in Fig. [51 3 of which are shown in Fig. [HI The correspond- 
ing identifications are listed on the last line of Table HI 
which are reasonable given the source frequency, and the 
frequencies of the resonances. The positions of the test 
sources were moved randomly, and it was always found 
that 5 source positions were needed to recover all of the 
dynamical modes accurately. If more than 5 were used, 
the recovered modes and eigenvalues did not change, as 
would be expected for a well-conditioned problem. The 
conclusion is that a simple interferometric measurement, 
at a single frequency, based on measuring the total power 
absorbed, has enabled the forms and responsivities of the 
dynamical modes to be found. 

Figure Uni shows the total power absorbed, in arbi- 
trary units, as a function of frequency for 11, z-polarized 
dipoles, placed 0.1 mm apart, each having ujo = 300 GHz, 
and F = 20 GHz. All of the dipoles except the central 
one had a polarizability of a/eo = 0.0005 mm"^, whereas 
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TABLE I. Fringe parameters of the 5 dipole structure de- 
scribed in the text. Two 2-polarized, 280 GHz, sources were 
placed in pairs at the positions x =1.0, 2.0, 3.0, 4.0, and 
5.0 mm, with y — 5.0 mm in all cases. The top block of data 
gives the amplitudes of the elements of the H matrix, and 
the middle block gives the corresponding phases. All values 
were recovered from the 'measured' fringes. The penultimate 
line gives the eigenvalues of the natural modes at the dipoles, 
and the bottom line gives the frequency (GHz) at which the 
corresponding mode is resonant. 
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FIG. 10. Total power absorbed (arbitrary units) as a function 
of frequency for 11 parallel dipoles, 0.1 mm apart, each hav- 
ing uuo ~ 300 GHz, and F = 20 GHz. All of the dipoles except 
the central one had a polarizability of a/eo = 0.0005 mm^, 
whereas the central one had (solid) a/eo = 0.001 mm"^, (dot- 
ted) a/eo = 0.003 mm^, and (dashed) a/eo = 0.005 mm^. 



the central one had (solid) a/eo = 0.001 mni'^, (dotted) 
a/eo = 0.003 mm^, and (dashed) a/eo = 0.005 mm^. The 
z-polarized source was located at x=0 and y—1 mm, and 
operated in the near-field mode. Even when the polar- 
izability is low, there is evidence of interaction, because 
the peak is at 310 GHz, not 300 GHz. As the polarizabil- 
ity of the central dipole is increased, the spectrum splits 
into strong features at 276 GHz and 323 GHz, indicating 
the excitation of two resonant modes. Figure [11] shows 
the form of the strongest mode at each of 276, 300, and 
323 GHz. At 276 GHz, the highly polarized dipole in- 
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FIG. 11. (Color online) Modal forms of the strongest ab- 
sorptive modes at (solid) 276 GHz, (dashed) 323 GHz, and 
(dotted) 300 GHz of the dashed spectrum shown in Fig. 1101 
All dipoles are either perfectly in phase or perfectly out of 
phase, and the lines are drawn to guide the eye only. 



duces a localized excitation, with the nearest neighbors 
180° out of phase with the central dipole. At 323 GHz, 
the mode is still localized, but now the nearest neighbors 
are in phase with the central, exciting dipole. There is no 
spectral feature at 300 GHz, and the dominant mode is 
less localized, with the nearest neighbors being station- 
ary, and the next nearest neighbors being out of phase. 
The excitation is influenced by the boundary conditions 
at the ends of the chain. The modal forms change sig- 
nificantly over the 280-320 GHz range, and can take on 
forms where the central dipole launches traveling, albeit 
rapidly decaying, waves. Beating can also be seen be- 
tween different spatial frequencies. 

Figure [T^ shows the amplitudes and phases of the con- 
tinuous visibilities when two sources are swept along a 
parallel line 1 mm from the chain. The sources were op- 
erated in near-field mode, and the plots correspond to 
frequencies of (solid) 276 GHz, (dashed) 323 GHz, and 
(dotted) 300 GHz. The fringes give a clear indication of 
the changing form of the localized mode. A set of visibil- 
ities were recorded with only 11 source positions, at the 
same distance and spanning the spatial range shown in 
Fig. [T21 It was found that at each frequency the forms of 
all of the dynamical modes and their responsivities could 
be recovered, regardless of whether the probes were op- 
erated in the near-field or far-field mode. As the total 
number of source positions was increased from 2, the re- 
covered modes were in error until 11 source positions were 
used, at which point all of the modal forms were found 
precisely. As the number of source positions increased 
above 11, no change was seen, showing that convergence 
could be reached, and recognized. 

FigurelT^shows fringe visibilities at 300 GHz, when the 
position of the highly polarizable dipole was shifted one 
place at a time: (solid) x =0.0 mm, (dashed) x =0.1 mm, 
(dashed) x =0.2 mm. The sources were swept along a 
line 1 mm from the chain of dipoles, and the visibilities 
referenced to the central position. The fringe is sensitive 
to the position of the localized mode. Moreover, all of 
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FIG. 12. (a) Amplitude and (b) phase of the recorded fringe 
visibility at frequencies of (solid) 276 GHz, (dashed) 323 GHz, 
and (dotted) 300 GHz. The sources were swept along a line, 
1 mm from the chain of dipoles, and the visibilities are refer- 
enced to the central position. 
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FIG. 13. (a) Amplitude and (b) phase of the recorded 
fringe visibility at 300 GHz. A highly polarizable dipole was 
placed at (solid) x =0.0 mm, (dashed) x =0.1 mm, (dashed) 
X =0.2 mm. The sources were swept along a line, 1 mm from 
the chain of dipoles, and the visibilities are referenced to the 
central position. 
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FIG. 14. (Color online) Modal forms of the strongest ab- 
sorptive modes at (solid) 270 GHz, (dashed) 330 GHz, and 
(dotted) 284 GHz of a system comprising 11 dipoles. All of 
the dipoles had a polarizability of a/eo = 0.0005 mm'^, except 
the ones at x =-0.2 mm and x =+0.2 mm (dipoles 4 and 8), 
which had a polarizability of a/eo ~ 0.005 mm"^. 
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FIG. 15. Total power absorbed (arbitrary units) as a sin- 
gle probe is scanned along a line parallel to a linear chain 
of 21 dipoles placed 0.2 mm apart, having luq =300 GHz, 
r =20 GHz, and a/ eg =0.5 mm'' The scan distance was 1 mm. 
The lines correspond to frequencies of (dotted) 42 GHz, 
(dashed) 205 GHz, and (solid) 500 GHz: the powers have 
been divided by 5x10*, 1x10^, and 1x10^ respectively be- 
fore plotting. 

the modal forms could agam be recovered with 11 source 
positions, allowing the position of the highly polarisable 
dipole to be located. 

More complicated linear systems were then considered. 
For example. Fig. [14] shows the modal forms of two in- 
teracting localized excitations. All of the dipoles had a 
polarizability of a/ep = 0.0005 mm^, except the ones at 
X =-0.2 mm and x =+0.2 mm (dipoles 4 and 8), which 
had a polarizability of a/eo — 0.005 mm^. The two highly 
polarized dipoles now oscillate either in phase or out of 
phase, with the interaction mediated by the background 
chain. Again, we were able to recover all of the modes 
and responsivities at any frequency using EAI. 

Finally consider a linear array of 21 dipoles, placed 
0.2 mm apart, having wq =300GHz, F =20 GHz, and 
a/eo =0.5 mm^. The spectrum, not shown, has high- 
Q resonances at 42 GHz and 205 GHz, each of which 
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FIG. 16. Amplitude of the continuous fringe visibility for 
an array of 21 dipoles placed 0.2 mm apart, ojq =300GHz, 
and F =20 GHz. (a) 42 GHz, (b) 205 GHz, (c) 510 GHz. 
The reference positions are (solid) mm, (dashed) -2.5, and 
(dotted) -5.0 mm in all cases. The ends of the chain are at 
±2 mm. 



has a single dominant eigenvalue. Figure [15] shows the 
power absorbed as a single probe was swept parallel to 
the chain, 1 mm from it, thereby producing a near-field 
scan. The two sweeps at (dotted) 42 GHz and (dashed) 
205 GHz, show strong evidence of end effects, whereas the 
sweep at (solid) 500 GHz, where no spectral features are 
present, has a form typical of a continuous nanowire or 
sheet of resistive material. Figure[T6]shows the amplitude 
of the fringe visibility at (top) 42 GHz, (middle) 205 GHz, 
(bottom) 500 GHz. The reference positions are (solid) 
mm, (dashed) -2.5, and (dotted) -5.0 mm in all cases. 
The ends of the chain are at ±2 mm. 

Using pairs of dipoles at the same scan distance, we re- 
covered the dynamical modes using EAI, and compared 
them with those calculated directly. A limited number 
of test-source positions, typically 20, were used, but no 
differences were seen when over 1000 were used. The 
behaviour at 42 GHz could always be recovered easily, 
and the behaviour at 205 GHz and 510 GHz could be 
recovered for linear scans <1 mm from the chain. At 
greater distances, recovery became uncertain, because of 
difficulty exciting dipoles near the ends. By always plac- 
ing 2 sources in line with the dipoles, but beyond the 
ends, the problem could be solved. The situation is sim- 
ilar to edge absorption in continuous films. Figure [T7] 
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FIG. 17. (Color online) Modal forms of the strongest absorp- 
tive modes at (dotted) 42 GHz, (dashed) 205 GHz, and (solid) 
500 GHz of a linear chain of 21 dipoles having placed 0.2 mm 
apart, having luq =300GHz, F =20 GHz, and a/eo =0.5 mm''. 
The lines are drawn to guide the eye only. 



shows the recovered modes in those cases where suffi- 
cient source positions were used to recover the modes ac- 
curately. The mode at (dotted) 42 GHz is a quasi-static 
solution, dominated by near-field coupling; where half of 
the dipoles are stationary, and the other half oscillate 
symmetrically, in perfect antiphase, about the centre of 
the chain. The fringe visibility, plotted in the top pane of 
Fig. [TBI shows directly that dipole interaction is strong 
enough to produce highly coherent extended behaviour 
along the whole length of the structure. Near-field inter- 
actions on linear chains are being considered for possible 
application as sub-wavelength waveguides. The mode at 
(dashed) 205 GHz is mediated by far-field retarded in- 
teraction. The dipoles at the ends are hghtly bound, 
because they only have neighbors on one side, and os- 
cillate with high amplitude; they form waves that travel 
inwards and interfere to produce a standing wave at the 
centre. The visibility plots shown in the middle pane of 
Fig. [16] reveal the existence of this large-scale collective 
behaviour. The bottom pane of Fig. [THj shows the vis- 
ibility at 500 GHz, where no obvious spectral features 
are present. The two largest eigenmodes are degenerate, 
each corresponding to a high amplitude end wave that 
travels towards the centre: only one of them is shown in 
Fig. [T7| (solid). Crucially, however, the effects at oppo- 
site ends are not coupled, as they are at 205 GHz. This 
behaviour is apparent in the visibility plots. Near the 
centre, the system behaves as an almost incoherent ab- 
sorber, whereas at the ends, highly coherent absorption 
is seen. The absorption process at one end is incoher- 
ent with respect to the absorption process at the other. 
Compare the bottom pane of Fig. [THj with the middle 
and top panes. The various forms seen are typical of 
thin-film planar absorbers, and therefore the linear chain 
is becoming large enough to take on the characteristics 
of a continuous absorber, and macroscopic multipolar ex- 
pansions can be used to describe behaviour (ll| . 

Many other simulations were performed during the 
course of the work. For example, we looked at arrange- 



ments where a single, lightly coupled dipole was placed 
next to a ring-like or linear chain, and found that we were 
able to separate out the behaviour of the lightly bound 
dipole from the collective behaviour of the chain. Also, 
we could identify the existence of 'grain-boundary-like' 
defects on long chains. 



IV. PRACTICAL CONSIDERATIONS 

Although this paper focuses on the theory of EAI, it 
beneficial to comment on its practical application. The 
central requirements are (i) the ability to phase-lock two 
sources, (ii) the ability to scan the sources within a near- 
field or far-field volume, (iii) the ability to measure the 
total power absorbed, or at least the part of it that con- 
stitutes the operation of the device. Because the key 
elements are so simple, the scheme can be used at any 
wavelength. Crucially, high-power sources can be used, 
perhaps at low temperatures, extending the method well 
beyond anything that can be achieved bymeasuring the 
correlations in thermally radiated fields |33l-l35j. 

It is certainly straightforward to assemble phased- 
locked sources for radio, infrared, and optical wave- 
lengths. In fact, we have already demonstrated that it 
is possible to measure the state of coherence of the re- 
ception patterns of multimode power detectors over the 
range 195-270 GHz using absorbed power only [s^. In 
that work, the frequencies of the two sources were offset 
slightly to produce a time varying output that contained 
the fringe information in real time. In some cases, a 
single-frequency source would be used, say in the case 
of a laser and power splitter, and then A/4 path- length 
switching would be needed to determine the quadrature 
components of the fringe. 

There are many ways in which the SUT could be illu- 
minated with different basis sets to span the dynamical 
modes of interest. For example, far- field probes in the 
form of optical fibres, or near-field probes in the form 
of coaxial AFM-like tips, could be scanned across or ro- 
tated around the SUT. In magnetic measurements, coils 
could be swept around and scanned across the surface 
of interest. In general, polarisation switching is needed, 
but this can be done in a straightforward manner. A 
crucial question is 'how many samples must be taken to 
gain access to the degrees of freedom of interest?' We 
have already identified efficient real-time techniques for 
adding in more and more sample points until convergence 
is reached, and the significant degrees of freedom identi- 
fied. EAI is, in general terms, quite similar to Aperture 
Synthesis Interferometry, used extensively in astronomy, 
and many of the data processing techniques developed 
for astronomical imaging, such as closure phase, could 
be incorporated here. 

The final requirement is a way of determining the to- 
tal power absorbed. In many applications, such as pho- 
tovoltaics and optical sensors, the output of the device 
is already a direct measure of the power being absorbed 
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through the degrees of freedom of interest. In these cases, 
one might not even be interested in the surface physics 
per se, but a detailed characterisation of the nature, num- 
ber, and responsivities of the electromagnetic modes to 
which the device is sensitive. Additionally, representa- 
tive models based on DDA can be used to understand 
how the electromagnetic modes relate to the dynamical 
behaviour of the absorbing surface. In some cases, such 
as patterned films or biomolecules on surfaces, a way of 
monitoring the total power absorbed is required. Our 
own preferred approach is to attach the SUT to a micro- 
machined SiN membrane, which also supports a thin-film 
temperature sensor. Such devices have been developed 
extensively for astronomical detectors, and when used at 
low temperatures can achieve extreme sensitivity j37l[38l |. 
Our detectors [soj can detect absorbed powers of much 
less than 1 fW with noise levels of 1x10"^^ WHz^^/^ 
making extreme thermometry of thin-film structures en- 
tirely possible. 

V. CONCLUSIONS 

We have shown that by illuminating an energy- 
absorbing surface or structure with a pair of phase-locked 
sources, it is possible to measure the state of coherence 
of the electromagnetic field to which the structure is sen- 
sitive. The tensor that describes the coherence can be 



diagonalized to find the amplitude, phase, and polar- 
ization patterns of the individual electromagnetic modes 
through which the structure can absorb power. In some 
cases, such as photovoltaics and detectors, this informa- 
tion is already sufficient to understand and characterize 
behaviour. It is also possible, however, to map the elec- 
tromagnetic modes back onto the surface, or object, to 
gain representations of the dynamical modes responsi- 
ble for absorption. These modes are intimately related 
to the existence of coherent phenomena, such as elastic 
waves and surface plasmons. In a continuous formula- 
tion of the method, not yet published, we have shown 
that it is possible to detect the presence of subsurface 
voids. The method can be implemented easily, and be- 
cause high-power sources are used, can far outperform 
any experiment that relies on measuring the correlations 
in themally radiated fields. More speculatively, the use 
of high power sources, together with a space-frequency 
formulation, may open the way to studying non-linear 
surface effects. 
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